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A comparative study of the magnetic remanence exhibited 
by tridimensional Josephson junction arrays in response to 
an excitation with an AC magnetic field is presented. The 
observed temperature behavior of the remanence curves for 
disordered arrays fabricated from three different materials 
{Nb, Y Ba2Cuz07-s and Lai,%z>Sro.\zCu04,-s) is found to 
follow the same universal law (based on the explicit tempera- 
ture expressions for the activation energy and the inductance- 
dominated contribution to the magnetization of the array 
within the framework of the phase-slip model) regardless of 
the origin of the superconducting electrodes of the junctions 
which form the array. 

PACS numbers: 74.50.+r; 74.25. Ha; 74.60. Jg 

As it has been recently found [0-^, tridimensional 
disordered Josephson junction arrays (JJAs) fabricated 
from either low-Tc (LTS) or high- Tc (HTS) superconduc- 
tors may, upon excitation by a magnetic field, exhibit 
a temperature-dependent magnetic remanence, Mji[T). 
Typically , the magnetized state occurs in a rather nar- 
row window of temperatures, the extent of which depends 
on the critical current, Ic{T), of the junctions. Besides, 
there is a threshold value for the magnetic field in order 
to drive the JJA to the state where flux is retained after 
suppression of the field Q . 
I In this Letter we present a comparative study of three 
different samples with a rather spectacular remanent be- 
havior and suggest a possible interpretation of the ob- 
served temperature dependence of the remanent magne- 
tization of both LTS and HTS tridimensional disordered 
JJAs. Our analysis shows that all the experimental data 
I can be rather well fitted using the explicit temperature 
expressions for the activation energy and the inductance- 
dominated contribution to the magnetization of the array 
within the so-called phase-slip model |^-^ . 

Three samples were prepared from selected mate- 
rial, respectively of Nb, Y Ba2Cuj,0'!^s (YBCO) and 
Lai,%^SrQ,i^CuOi-s (LSCO). All three exhibit the pre- 
dicted remanence and other characteristic features of 
Josephson arrays. Fabrication routes as well as the ex- 



perimental routines employed for the magnetic measure- 
ments are described elsewhere In short, the corre- 
sponding (e.g., niobium) powder was separated accord- 
ing to grain size (using a set of special sieves, with mesh 
gauges ranging from 38 to 44/xm), then uniaxially pressed 
in a mold to form a cylindrical pellet of 2.5mm radius 
by 2.0mm height. This pellet is a tridimensional dis- 
ordered JJA in which the junctions are weakly-coupled 
grains, i.e., weak- links formed by a sandwich between 
(Nb) grains and a (iV&-oxide) layer originally present on 
the grain surface. 

The measurements were made using a Quantum Design 
MPMS-5T SQUID magnetometer featuring the regular 
DC extraction magnetometer and an AC-susceptibility 
module. The remanence was obtained measuring the 
sample magnetization after application and removal of 
a train of sinusoidal pulses. Using the field scan routine 
we measured the remanent magnetization as a function 
of the excitation field. For an ordinary superconductor 
of any kind, from a single crystal to a totally disordered 
granular sample, the only possibility of a remanence af- 
ter the application of the AC field would be a residual 
magnetization due to flux eventually pinned inside the 
specimen. This contribution, however, is expected to 
be small and practically independent of the excitation 
field. We have verified the above characteristics measur- 
ing Mfi{h,T) for a variety of samples. In particular, the 
powder used to fabricate our arrays have the typical re- 
sponse of ordinary superconductors, so that the effects 
described below are entirely due to the formation of the 
3D- JJA. The experimental results for all three samples 
(along with the model fits, see below) are summarized in 
Fig.l which suggests that the observed behavior seems 
to follow a universal temperature pattern, irrespective of 
the type of superconductor of which the array is made. 
Let us turn to a possible interpretation of the obtained 
results. 

Since the observed remanent magnetization (RM) in 
our samples (JJAs) appears only below the so-called 
phase-locking temperature Tj (which marks the estab- 
lishment of phase coherence between the adjacent grains 
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in the array and always lies below a single grain super- 
conducting temperature Tc), it is quite reasonable to as- 
sume that origin of RM is related to thermal fluctua- 
tions of the phases of the superconducting order parame- 
ters across an array of Josephson junctions (the so-called 
phase-slip mechanism ||-^)- In the present approach 
we consider the sample as a single plaquette with four 
Josephson junctions (JJs), each of which is treated via 
an effective single junction approximation. Within this 
approximation, the phase-slip scenario yields then 

AMr{T) = M{T) - Mr = Mo{T)I^^Yf{T)/2] - Mr 

(1) 

for the observed remanent magnetization. Here, Mq(T) is 
an inductance-induced contribution to the magnetization 
of the array (see below), 7(T) = U{T)/kBT is the nor- 
malized barrier height for thermal phase slippage, /o(a;) 
is the modified Bessel function, and Mr = M{Tj) is 
a residual temperature-independent contribution (notice 
that, according to Eq.(l), AMr{Tj) = 0). 

For temperatures below Tj (where the main events 
take place, see Fig.l), the Bessel function can be ap- 
proximated as loix) ~ l\J1-KX leading to a simphfied 
version of Eq.(l): 

M{T) = 2^A/o(r)[;7(r)/fcsr] exp[-;7(r)/fcsr] (2) 

Figure 1 shows the temperature dependence (in reduced 
units, T = T jTj) of the normalized remanent magnetiza- 
tion mr(r) = I\Mr{T)I 1\Mr{Tp), where Tp is the peak 
temperature and AMij(T) is defined via Eqs. (1) and 
(2). The data for YBCO- and Nb-based J J As are found 
to be well fitted with the following explicit expression for 
the array magnetization: 



M{t) = A(\ - t'^fl'^ exp[-a(l - t% 



(3) 



where t — T/Tc- The best fits through all the data points 
(shown in Fig.l by solid and dotted lines for YBCO- and 
Nb-based JJAs, respectively) using Eq.(3) and the known 
critical parameters 

YBCO : Tc = 90K, T, = 82K, Tp = OMT,; (4) 



Nb 



Tc^9.1K, Tj^8.2K, Tr = 0.92Tj; (5) 



yield the following estimates of the model parameters: 
avBCO = 7 and aNb = 9. 

At the same time, the data for the LSCO sample 
(which appears to have links weaker than the other two 
samples) are found to be better fitted using the following 
expression; 



M{t) = B{1 - ^2)5/2 exp[-/3(l - ^2)3/2]^ 



(6) 



LSCO : Tc = 36.5K, Tj = ig.STif, Tr = 0.7Tj; 

(7) 

yielding hsco = 2. 

To understand the observed behavior of the rema- 
nent magnetization, we need to specify the tempera- 
ture dependencies of the activation energy U{T) and the 
inductance-dominated contribution AIq(T) to the mag- 
netization of the array. Starting with the YBCO- and 
Nb-based arrays, it is reasonable to assume that [pO 11 
U{T) = $o-fc(7')/27r and Mq{T) = LIc{T)/hqS, where 
Ic{T) is an average value of the critical current, L is 
an average inductance of the Josephson network, S is an 
effective (in general, temperature-dependent, see below) 
projected area of the contact, $0 is the flux quantum, 
and /io is the vacuum permeability. In turn, the tem- 
perature dependence of the critical current is dictated by 
the corresponding dependence of the London penetration 
depth, namely M: 



Ic{T)^IcmXLiO)/XL{T)]' 



where 



AL(r) = Ai(0)[l 



(T/Tc)^]-^/' 



(8) 



(9) 



Finally, to arrive at the fitting expression given by 
Eq.(3), we have to assume that the projected area S 
is also temperature dependent (which is not unusual), 
viz. S{T) = iid[T)l with d{T) and I being the thick- 
ness and the length of a SIS-type sandwich, respectively 
(d(T) = 2\l{T)+^, where Al(T) is the London penetra- 
tion depth and ^ is the thickness of an insulating layer; 
in ceramics / plays the role of an average grain size r^; 

typically, Z»Al(T)» 0- 

Turning now to the LSCO-based array, we notice that 
in this case the temperature dependence of the activation 
energy U(T) is governed by the well-known weak- links 
mediated expression g U{T) = C/(0)[1 - {T/Tcff''^, 
while the critical current reads Ic{T) = /c(0)[l — 
(T/Tcf] (Cf. Eq.(8)). Besides, in this particular case 
the projected area is temperature-independent, S = ttP. 

The above considerations bring about the following re- 
lationships between the fitting and the model parametres: 



Llc{0)a 
/ioAL(0)r 



B = 



2LIciO)l3 



(10) 



with 



with a = <S>oIc{0)/2TTkBTc and /3 = U{0)/kBTc. 

In conclusion, to check the self-consistency of the 
model, we recall that within our scenario, the normalized 
remanent magnetization mr{t) disappears for T > Tj 
(see Fig.l) where Tj is the phase- locking temperature of 
the Josephson network. At the same time, this temper- 
ature is usually defined via the equation U{Tj) = ksTj 
which results in the following two expressions (both valid 
near Tc) relating the two critical temperatures {Tc and 
Tj): 
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and 



{l-tjf'^ = 



tj_ 



(11) 



(12) 



for YBCO and Nb, and for LSCO based arrays, respec- 
tively. Using the above-mentioned values of the critical 
temperatures for the three samples and the experimen- 
tally found estimates for a and f3 (sec above), Eqs.(ll) 
and (12) bring about the following reasonable estimates 
of the phase-locking temperatures (Cf. Eqs.(4), (5) and 
(7)) for the three arrays: Tf^C'O = (7/8)T^^<=''' , 

rpLSCO ^ {l/2)T^SCO r^Nb ^ (12/13)^^6. 

In summary, by employing the so-called phase-slip 
mechanism in each clement of a plaqucttc with four ef- 
fective Josephson junctions, we have consistently derived 
expressions which lead to a reasonable description of the 
temperature behavior of magnetic remanencc of disor- 
dered Josephson arrays. The values obtained for the two 
experimentally accessible parameters, the locking tem- 
perature, Tj, and the superconducting critical tempera- 
ture, Tc, are in good agreement with the experimental 
data. The employed approach indicates that the temper- 
ature dependence of the magnetic remanencc is universal, 
regardless of the origin of the superconducting electrodes 
of the junctions which form the array. 
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FIG. 1. Temperature dependence of the normalized rema- 
nent magnetization mr(T), showing the experimental data 
for three different samples and the corresponding fittings (see 
text). 
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